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1 Introduction 



In this article, we study the global well-posedness of the semi-relativistic Schrodinger-Poisson 
system on a finite domain. This system is relevant to the description of many-body semi- 
relativistic quantum particles in the mean-field limit (for instance, in heated plasma), when 
the particles move with extremely high velocities. Consider semi-relativistic quantum parti- 
cles confined in domain C which is an open, finite volume set with a boundary. The 
particles interact by the electrostatic field they collectively generate. In the mean-field limit, 
the density matrix that describes the mixed state of the system satisfies the Hartree-von 
Neumann equation 

'idtp{t) = [Hv,p{t)], xeQ, t>0 ^^^^ 
-AV' = n(t,x), n{t,x) = p{t,x,x), p(0) = po 

satisfying Dirichlet boundary conditions, p{t, x,y) = if x or y E dQ, for t > 0. The 
Hamiltonian is given by 

Hv:=T^. + V{t,x) (1.2) 



where the relativistic kinetic energy operator '■= \^—A + m? — m is defined via the 
spectral calculus. Here, A denotes the Dirichlet Laplacian on L^(f2), and m > is the 
particle mass; see [H |2] for a derivation of this system of equations in the non-relativistic 
case. Since p{t) is a positive, self-adjoint trace-class operator acting on L^(f2), its kernel can, 
for every t G M+, be decomposed with respect to an orthonormal basis of L'^{VL). The kernel 
of the initial data po can be represented in the form 

Po{.x, y) = ^ \kilJk{,x)i)k{y) (1.3) 

where {il)k}k&i denotes an orthonormal basis of L^(f2), with 'ipk\dn = for all k E N, and 
coefficients 

X:={h}kmei\ Xk>0, 5^Afc = l. (1.4) 

As shown below, there exists a one-parameter family of complete orthonormal bases of L'^{Q), 
{'0jfc(t)}fcgN5 with 'ipk(t)\sn = for all k E N, and for t G M+, such that the kernel of the 
solution p{t) to (11. ip can be represented as 

p{t, x,y) = Xki^k{t, x)tlJkit, y). (1.5) 



Notably, the coefficients A are independent of t, and thus the same as those in po- Substituting 
(11. 5p in (II. ip . the one-parameter family of orthonormal vectors {'ipkit)}km is seen to satisfy 
the semi-relativistic Schrodinger-Poisson system 

= Tmi^k + ViJk, keN (1.6) 
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-AVm=n[% vl/:={^,}-^, (1.7) 

oo 

n[^{x,t)] = Y,h\i^k\', (1.8) 

k=l 

with initial data {"0^(0)}^;^. The potential function ^[^1^] solves the Poisson equation fll.7p . 
On both V[^] and ipk{t), for all ^ G N, we impose Dirichlet boundary conditions 

ipk{t,x) , V{x,t) = 0, t>0, e dQ. (1.9) 

As we show in Lemma El below, solutions of fll.6p - fll.8l) preserve the orthonormality of 
{V^fc(^)}fceN- 

The state space for the Schrodinger-Poisson system is given by 

1 

^ '■= {(*,A) I * = {'4'k}'k=i C H^i^)nH^{^l) is a complete orthonormal system in L^(fi), 

oo „ 

X = {Xk}k=i^£\ Xk>0, keN, 5^Afc / |V^fcp(ix<oo}. 

k=i -l^ 

For fixed A G > 0, and for sequences of square integrable functions $ := {(t)k}T=i 

^ := {V^fcjfcLii we define the inner product 

oo 

■>'^) cl{Q.) ■= ^Xk{(j)k,tpk)L^n), 

k=l 

which induces the norm 

oo 

ll*ll£|_(n) := Afc||0fc||^2(n))2, 

k=l 

and we introduce the corresponding Hilbert space 

4(0) := {$ = 1 I 0, e L'{n), V A; G N, mclm < oo}. 

Our main result is as follows. 

Theorem 1. For every initial state (\&(x,0), A) G £, there is a unique mild solution '${x,t), 
t G [0,oo), of U.6\) - in~^) with (\E'(a;, t), A) G C, which is also a unique strong global solution 
m £|(fi). 

Establishing the global well-posedness of the Schrodinger-Poisson system plays a crucial 
role in proving the existence and nonlinear stability of stationary states, i.e. the nonlinear 
bound states of the Schrodinger-Poisson system, which was done in the nonrelativistic case 
in [U [6]. The problem in one dimension was treated in [8]. The semiclassical limit of the 
Schrodinger-Poisson system with the relativistic kinetic energy was studied in the recent 
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article [T]. Global well-posedness for a single semi-relativistic Hartree equation in was 
established in [5]. In the present work, we deal with the infinite system of equations in a 
finite volume set with Dirichlet boundary conditions, and, as distinct from [5], we do not use 
the regularization of the Poisson equation. Moreover, both the results of [5] and Theorem [1] 
above do not rely on Strichartz type estimates. 



2 Proof of global well-posedness 

We make a fixed choice of A = {Xk}'^^i G with > and ^A^ = 1, denoting the 
sequence of coefficients determined by the initial data po of the Hartree- von Neumann equa- 
tion fll.ip via fll.Sp . for t = 0. We note that we require all A^ > to be positive for the 
subsequent analysis. This does not lead to any loss of generality since by density arguments, 
any po (and likewise p{t)) can be approximated arbitrarily well by an expansion of the form 
fOj) . respectively ([13]), with A^ > 0. 

We introduce inner products (') O-^i/a^^^^ and {■,-)nl{n) which induce the generalized in- 
homogenous Sobolev norms 

oo oo 
~ k=l k=l 

and define the corresponding Hilbert spaces 
and J 

respectively. We also introduce the generalized homogenous Sobolev norms 



•= ^k\\\p\^k\\l2(^n))^ and := (J^ Afc|| V(/)fe||i2(f^) 

k=l k=l 



2. 



Here, \p\ stands for the operator y— A? and has the meaning of the relativistic kinetic energy 
of a particle with zero mass. We note the following equivalence of norms. 

Lemma 2. For $ G 1-L]/^{^1), the norms ll^ll-^va^-f^-) and II 1 1 ^1/2 ^-f^^ are equivalent. // $ G 
'H\{Q), then W^WnKn) equivalent to ||"^'||-^i(n) ■ 
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Proof. Clearly 

oo oo 
k=l k=l 

We will make use of the Poincare inequality 

[ \V<Pk\^dx >Cp I \(t)k?dx (2.1) 

with the constant Cp > dependent upon the domain Q with Dirichlet boundary conditions. 
Thus 

|2 \ l|2 



which enables us to estimate 



k=l 



V k=i - 

Let us compare the remaining two norms. Clearly, 

oo 
k=l 

On the other hand, by means of the Poincare inequality fl2.ll) . 

oo 



k=l 



< 



II 1 

k=i 

□ 

Let = {'0m}m=i be a wave function and the relativistic kinetic energy operator acts on 
it Tm"^ = (v— A + ?7? — m)ilj componentwise. We have the following two lemmas. 

Lemma 3. The domain of the kinetic energy operator is given by D{Tm) = 'H\{Q) C C'^[Q). 
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Proof. Let * G n\{n). Then 



oo oo oo 

A^IIV'mllHHQ) = ^rniWi^mWh^n) + II || i2(j^) } > J] ^^117/^^1112(^2), 

m=l m=l m=l 



and also, 11^11^2 (j^) < oo. We estimate 



Il^m'0fc||i2(f2) = ((-A + m^)?/^fc,^fc)L2(o) + "^^||V'fe|li2(f2) -2m{V^ATm?'ilJk,ipk)L^n) < 

< i|VV^fc||i2(f2) + 2771211^^^112 2^^^ < c(m) II V^fc II ^i(f^), 
where c(m) is a mass dependent constant. Hence 



k=l k=l 



l^'n^llrKn) = ^kWTmiJkWhf^n) < cM ^ Afc||7^fc||^i(f2) < oo. 

□ 



Lemma 4. T/ie operator generates the group e ^^^^ t G M, o/ unitary operators on 

Clin). 



Proof. For a, /3 G /^^(fi) we compute the inner product 

oo oo 
k=l k=l 

□ 



We rewrite the Schrodinger-Poisson system for x E fl into the form 

*t = -iTm"^ + t)], Tw/iere F[^] := rV[^]^, (2.2) 

-AF[^] = n[^], w;/iere V|an = 0, 

oo 

fc=i 

and prove the following auxiliary result. 

Lemma 5. The map defined in Ii2.2\) F : 'H\{Q) — t- 'H\{Q) is locally Lipschitz continuous. 
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Proof. Let ^, $ G H\{n) with * = {^4^=1, $ = {04^1 and t e [0,r]. Then, 

iiF[vi>]-F[$]|i„^(^) = \\t~'vm^-z~'v[m\Him = \\vmi<ii~^)+ivm-v[^mk^^^n)- 

This can be easily estimated above by means of Lemma [2] by 

c\\vmi^ - mnim + c|l(v^[^] - vmMnim^ 

which equals 

oo oo 

C{Y.^>'\\^(nmi^k-A))\\l.^n)y' + C{J2^^^^^ (2-3) 

k=l k=l 

Here, C denotes a finite, positive, universal constant. Clearly, we have 

<2||(Vl^[vl>])(^,-0,)|| 
By means of the Schwarz inequality this can be bounded above by 

C\\^Vm\\lHn)\\^k - M'mn) + 2\\Vm\\Uin)\\'^i^k - <l>k)\\h^ny 
By applying the Sobolev embedding theorems to these expressions, we arrive at 

To estimate the remaining term in (12. 3p . we use 

The Schwarz inequality yields 

2\\V{Vm - V[mlHn)m\lHn) + 2||(^[*] - V[mUin)\\'^MUnr 

Applying the Sobolev embedding theorem along with the Holder inequality to these expres- 
sions, we find 



C\\A{V[^] - F[$])||i.(^)||0fc||ia(^) + C\\A{Vm - V[mhin)mk\\h 



From the Sobolev inequality used in the first of the two terms above we deduce the upper 
bound 

c\\vm-vm\i,^^^\\v<p4 

L2(n)- 

Therefore, for the norm of the difference ||F['0] — F[<l>]||-^i(-f^) we have the estimate from above 

as 

oo oo 
k=l k=l 



which obviously equals to 

C\\VmHHn)\\'^ - nnlin) + C\\Vm - VmUnHnminim- 

Let us apply the Poincare and the Schwarz inequalities to estimate the Sobolev norm of the 
potential function as 

\\Vm\HHn) < C\\AV\\mn) = C\\nm\\LHn)- 

Hence, our goal is to estimate the appropriate norm of the particle concentration. From the 
Schwarz inequality, 



\n 



k,l=l k=l 

and using the Holder inequality along with the Sobolev inequality. 
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Hence, we arrive at the estimates for the particle concentration and the norms on the po- 
tential function, 

iin[^^]iu.(^) < cm\^^y mmnnn) < cm\^n) 

with II ■ 11-^1 and || ■ ||-Hi(n) equivalent via Lemma |2l Evidently, 

W := V[^] - 

satisfies the Poisson equation, 

-AW^ = n[^] -n[$], W\9n = 0, 

and Dirichlet boundary conditions. Applying the Poincare inequality along with the Schwarz 
inequality, we arrive at 

<C||Al^||i.(^), 

such that 

\\W\\H^n)<C\\n[^]-nmL^n). 
We will use the trivial inequality 

oo 

\n[^] - n[$]| < 5^ AfcdV^fcl + |0fc|)|^fc - 0fc|. 

k=l 



The Schwarz inequality applied twice yields 



oo 



< (J^ Afc|||V'fc| + |0fc|||L4(n)||V'fc-0fc||L4(n))^ < h{\\tpk\\L^n) + \\<Pk\\\LHn))Hk- (l)k\\LHn)y , 

k=l k=l 

and using it again gives 

oo oo 

^hiWipkh^in) + \\(t)k\\\L-i(n)f^K\\ips - 0.||i4(f^). 

k=l s=l 

Applying the Holder and Sobolev inequalities, we arrive at 

oo oo 

+ IIV0.II )^A,||V^. -V0.il 

k=l s=l 

This quantity can be easily estimated above by 

(oo oo \ oo 

k=l 1=1 J s=l 

which clearly equals to 

C(ll^ll?.i(n) + Il<^ll?.i(n))ll^-'^II?.^(C)- 

Therefore, 

||n[^] - n[$]|U2(n) < C(||^'||^^(f,) + ||$||Ki(n)) II* " '^'llKi(n) 

and 

\\y[^] - yMlHHu) < C(||^'||«^(f,) + 11$ II ^^(f,)) II ^' - nn\in)- 
Collecting the estimates above, we arrive at 

||F[vP] - F[$]||^.(^) < Cm\\^n) + m'n^^in)^ " nniin), 

which completes the proof of the lemma. □ 

From standard arguments (see for instance Theorem 1.7 of |7]) thus follows that the above 
Schrodinger-Poisson system admits a unique mild solution \Ef in 'H\{Q) on a time interval 
[0,T), for some T > 0, satisfying the integral equation 

^(t) = e-*^'"**(0) + e"*^'"* /" e'^^'F[^{s)]ds (2.4) 

Jo 



Ill nli^). Moreover, 

limi^T||^(Olk^(f^) = 

if T is finite. We also note that is a unique strong solution in C\{Q). We shall next prove 
that this solution is in fact global in time. First we prove the following lemma. 

Lemma 6. Suppose for the unique mild solution (12.41) of the Schrddinger-Poisson system 
U.6\) - I[l~^) that {tpklx, 0)}^]^ att = forms a complete orthonormal system in L'^{Q). Then, 
for any t G [0,T), the set {i(jk{x,t)}'^-^ remains a complete orthonormal system in L^(f2). 
Moreover, the C\{fl)-norm is preserved, \\'^{x,t)\\^2^(^Q-j = \\'^{x,0)\\^2^(^q-j, t G [0,r). 

Proof. Given the solution "^{t) of the Schrodinger-Poisson system on [0,T), we obtain the 
time-dependent one-particle Hamiltonian 

Hv^it) = T^ + V<,{t,x) 

where the potential V* solves — AV*(t,x) = n[\E'(t)] with Dirichlet boundary conditions, see 
(11. 2p . Accordingly, the components of solve the linear, non- autonomous Schrodinger 
equation idti/jkityx) = Hv^(t)'iljk{t,x), for G N, on the time interval [0,T). We thus have, 
for t G [0,T), 

M^^t) = (e-^'i'o^^*M^>fc)(x,0), keN, (2.5) 

and therefore 

{M^,t),Mx,t))L^n) = (e-^i'o^^*W'^>fc(x,0),e-*i'o^v,Mdr^^(3,^0))i2(o) = 

= {tpkix, 0),iJi{x, 0))L'2(n) = Sk,i, k,l eN, 
where 6k,i stands for the Kronecker symbol. Obviously, for A; G N, 

\\Mx,t)\\l2(^^) = \\Mx,o)\\l2(^n), 

such that for t G [0, T), the Cl{n) -norm is conserved, 

oo oo 

\\'^ix,mclin) = EAfc||^fc(x,t)||i2(f,))5 = (5^Afc||V'fc(x,0)||i.(^))5 = \\^{x,0)\\ci^n)- 

k=l k=l 

Let us consider an arbitrary function f{x) G L^{Q). Clearly, we have the expansion 

oo 

/(^) = ^ifiy)^'^k{y,0))L2(n)iJk{x,0) 

k=l 
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and similarly 

oo 
k=l 

Thus, by means of fl2.5p we arrive at the expansion 

oo 

fix) = ^{f{y),iJk{y,t))L^n)i^k{x,t) 

k=l 

ioTte[0,T). □ 

Furthermore, we have conservation of energy for solutions to the Schrodinger-Poisson 
system in the following sense. 

Lemma 7. For the unique mild solution i \2A\\ of the Schrodinger-Poisson system ( fi.gjj - pT^) 

and for any value of time t G [0, T) we have the identity 

+ l\\VV[^{x,t)]\\hin) = ll*(^'0)lliv.(^) + ^l|Vy[*(x,0)]||i.(^). 
Proof. Complex conjugation of the Schrodinger-Poisson system (I1.6P yields 

-i^ = r^^fc + y[V^]4, keN. (2.6) 

Adding the A;-th equation of the original system (11. 6p multiplied by , and the A;-th 

at 

dijj 

equation in (12. 6p multiplied by ^6 obtain 

§:\\TUk\\l2^n) + / V[ij]^\^k\'dx = 0, ke N. 
Thus, multiplying by A^, and summing over /c, we find 

^ll*(^,i)ll^V2(^) + _^^[*(a:,^)]^'^[^(3^,i)]'i^ = 0. (2.7) 
One can easily verify the identity 

||| Vy[vl/(x, f)]||i.(^) = 2 ^ y[vl/(x, t)]|n[vl/(x, t)\dx, 

which we substitute in (12. 7p to complete the proof of the lemma. □ 
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With the auxihary statements proven above at our disposal, we may now prove our main 
result, Theorem [U 



Proof of Theorem U\ The proof follows from the blow-up alternative and conservation laws. 
It follows from Lemma [7]that ||^(^)|Ll/2^o^ is bounded from above uniformly in time. 

We need to bound W'^it) ||-^i(n)- We recall the mild solution of the Schrodinger-Poisson system 
(OD-dlll), given by 

^(t) = e-^^'"*^(0) +e-*^'"* /" e^^'""F[^(s)](is, (2.8) 

Jo 

which implies 

||vl/(t)||^,(n)<||vl>(0)||^i(n)+ / ||F[^(s)]||^i(f,). 

-Jo 

From Lemma [21 we have 



\ 1/2 

<c[Y,h\\ViVmif;, 



^fcJllL2(f7) 

,fe=l 



Now, 

< l|VF[v&]||i6(^)||V^fci|i3(f,) + \\Vm\\Uin)\\'^i'k\\ L2(n) 
where we have used Holder's inequality in the second line and the Sobolev inequality 

II/II^^(^^<cii/i|h.(.) 

in the last line. To evaluate || Vy[^]||L6(Q), recall that AV[^] = Applying Holder's 

and Sobolev inequalities, we get 

oo oo 

k,l=l k,l=l 
oo oo oo 

m i/l/2(Q)) 
k,l=l k=l 1=1 

<CII\Ef|P-, ll^'ll^,/ 
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We now estimate || V^[^']||L°°(n)- The Sobolev inequality implies 

mmUn) < C\\\p\-'/Mn\lHny 

We claim that || ||L2(n) is controlled by ll^ll-^i/z^-f^-)- 



< C\mlHn)\\\p\~Mmm^Hn) < C\mm/2^n)\\\p\~'^HmLHn), 

where we have used Holder's inequality in the first line, and the Sobolev inequality in the 
second line. It follows that 

\\\pr/MmLHn)<cmi^/.^^y 



and hence 



Combining the above estimates yields 

l|i^[^]lki(.)<^^ll*ll^./.(,)ll*llHi(.)- 

This implies 

mm 



where Cq is a constant proportional to the initial energy ||^(0)||^. + |ll^^[^(0)]|li2cm- 



By Gronwall's lemma, 



<^^ie^^*, t>0. 



By the blow-up alternative, this implies that the Schrodinger-Poisson system is globally 
well-posed in 'H\{Q). □ 
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